Relations between the Gribov-horizon and 
center- vortex confinement scenarios 1 



Abstract. We review numerical evidence on connections between the center-vortex and Gribov- 
horizon confinement scenarios. 



Jeff Greensite*, Stefan Olejnfk 1 ^ and Daniel Zwanziger** 

* Physics and Astronomy Dept., San Francisco State University, San Francisco, CA 94117, USA 
^Institute of Physics, Slovak Academy of Sciences, SK-845 11 Bratislava, Slovakia 
** Physics Department, New York University, New York, NY 10003, USA 

> 

o 
co 

So oft in theologic wars, 
The disputants, I ween, 
Rail on in utter ignorance 

(NJ \ Of what each other mean, 

CO ■ And prate about an Elephant 

< ^ \ Not one of them has seen! 

John Godfrey Saxe (1816-1887) [1] 

o 



INTRODUCTION 



Once upon a time, several blind men from a Hindustan village touched different parts of 
an elephant's body and, judging from their perceptions, quarreled about what the whole 
creature might look like. In a similar way, we explore color confinement by isolating 
different aspects of the phenomenon, and build our limited experience into various model 
schemes. The moral from the ancient parable teaches us that it is vitally important to 
unify all different views before we can appreciate the whole beauty of the beast. 

The aspirations of this contribution are more modest. I will concentrate on common 
points of two seemingly unrelated pictures of color confinement. The former assumes 
that confinement arises due to the condensation of a particular type of topological 
excitations, so called center vortices (for a review see [2]), in the QCD vacuum; while 
the latter is the Gribov-horizon scenario in Coulomb gauge, which has been advocated 
by Gribov [3] and Zwanziger [4]. 

I will first briefly introduce the idea of the Gribov-horizon scenario, then formulate 
a simple criterion of confinement for static color charges in Coulomb gauge, discuss 
how the fulfillment of this criterion depends on presence/absence of center vortices, 
and finally present results for the Coulomb energy of a pair of static charges. I will 
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present a subset of our numerical results, details, as well as some analytic insights on 
the connections between center vortices and the Gribov horizon, can be found in recent 
publications [5, 6, 7]. 

CONFINEMENT SCENARIO IN COULOMB GAUGE 

In Coulomb gauge, the Hamiltonian of QCD has the following form [8]: 

H = H g i ue + H cou i, 
H 8 iue = ^A^-^EV-EV-^ + B^Bj, 

H coul = l -Jd\d 2 yf- l l 2 P a {x)fK ab {x,y-A)p b {y)f- ll \ 

P a = p a matt er-gf abc ± b -V. 
A prominent role in this expression is played by the Faddeev-Popov operator 

M(A) = -V-^(A), where @f c {A) = di8 ac + f abc A b (x), (5) 
which enters both the interaction kernel K and the Jacobian factor ^ 

K ab (x,y-A) = [M- 1 (-V 2 )M- ! ]^, J = det [M{A)\ . (6) 

It is well-known since the seminal paper of Gribov [3], that Coulomb-gauge fixing in a 
non-abelian theory is a difficult problem. The transversality condition V • A = does not 
fix the gauge completely. Gribov [3] suggested to restrict integration over configurations 
to the so-called Gribov region (GR), defined as the subspace of transverse gauge fields 
for which the Faddeev-Popov operator is positive, and which are therefore local minima 
of the functional 

I[A,g] = Jdx[8A a (x)} 2 , where *A i = g- 1 A i g + g- 1 d i g. (7) 

The boundary of this region is called the Gribov horizon. 

However, there exist gauge orbits which intersect the Gribov region more than once; 
the next step [9] then is to restrict fields to the fundamental modular region (FMR), the 
set of absolute minima of the functional (7). Both the GR and the FMR are bounded in 
every direction and convex. 

The essence of the Gribov-horizon confinement scenario can be phrased in a simple 
way: The dimension of gauge-field configuration space is huge, so it is reasonable to 
expect that most configurations are located close to its boundary (the horizon; in a 
similar way, the volume measure r d ~ l dr of a J-dimensional sphere is peaked at its 
surface). The interaction kernel K, which determines the interaction energy of static 
color sources, contains the inverse of the Faddeev-Popov operator, which is strictly zero 
on the horizon and near- zero close to the horizon. A high density of configurations near 
the horizon can thus lead to a strong enhancement of the Coulomb interaction energy, 
and hopefully cause color confinement. 
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A CONFINEMENT CONDITION IN TERMS OF 
FADDEEV-POPOV EIGENSTATES 



Let us consider a single static color charge, which can be written in Coulomb gauge as 

^[A;x] = xtr a (x)^ Q [Ai (8) 

where a is the color index for a point charge in color group representation r, and ¥0 is 
the Coulomb-gauge ground state. The excitation energy of this state, above the ground 
state, is given by 



Sr = (yg|yg) ~ W*"* W> (9) 

and is proportional to 

S =j^— [ (K aa (x,x;A)). (10) 

The quantity £ r is the color Coulomb self-energy of unscreened color charge and is 
expected to be both ultraviolet and infrared divergent in a confining theory. The UV 
divergence can be regulated by a lattice cut-off, however, the quantity must still be 
divergent at infinite volume, even after lattice regularization, due to IR effects. 

On a lattice, one can express (10) through eigenstates of the Faddeev-Popov operator 

YM a H n)h =K$ )a (id 

b,y 

simply as (assuming that M is invertible, i.e. excluding its zero modes) 

'=^l(^y wherc F n = £^ )a (-V 2 U4 n)a \ (12) 

-^3 n \ Arc / a,xy 

and V3 is the lattice 3-volume. In SU(2) lattice gauge theory, the link variables can be 
expressed as 

U^(x) = b^(x) + io c a c ^x), & M (x) 2 + J>£(x) 2 = l, (13) 

c 

and the lattice Faddeev-Popov operator 

< = d ab ^{d xy [b k (x)+b k (x-k)}-8 ' xy _M x )- 8 y, x -My)} 
k K } 

k k 

(x,y denote lattice sites at fixed time) is a 3V3 x 3V3 sparse matrix with 3V3 linearly 
independent eigenstates. If we denote N(X : X + AX) the number of eigenvalues in the 
range between X and X +AX, we can introduce, on a large lattice, the density of states 
p(X) =N(X, X +AX)/(3V 3 AX). Then, as V 3 -> 00, 

f hrnax rj 2 

£ = l j T (p(X)F(X)). (15) 



A (necessary) condition for confinement can now be formulated: The excitation energy 
$ r of a static, unscreened color charge is divergent if, at infinite volume, 



CENTER VORTICES AND THE CONFINEMENT CONDITION 

It is interesting to check whether the above condition is fulfilled in various ensembles of 
lattice configurations. First, at zero-th order in the gauge coupling, the Faddeev-Popov 
operator is simply a lattice Laplacian and its eigenstates are just plane waves. One can 
easily verify that, in this case, 

p(A) = ^v / I, F(X) = X, <g = ^\[l~^^ (17) 

the excitation energy is IR finite and the confinement criterion is not met. One needs 
some mechanism of enhancement of p(A) and F{X) in the region of small X values. 
I will argue that such an enhancement exists in full lattice configurations, and is provided 
by center vortices. 

Center vortices are identified by fixing to an adjoint gauge, and then projecting link 
variables to the subgroup of SU(/V) [10]. The excitations of the projected theory are 
known as P- vortices. In the direct maximal center gauge (DMCG) in SU(2) [11] one 
fixes to the maximum of 

^MCG = £|iTr[C/ M W]| 2 , (18) 

x,n 

and center projects by 

U^x) — Z M (x) = sign TrfC/^x)]. (19) 

A lot of evidence has been accumulated that center vortices alone reproduce much of 
confinement physics, for a review see [2]. 

We have determined the density of states p(A) and the mean value of the lattice 
Laplacian F(X) in three ensembles of lattice configurations: 

1. full configurations, {U^(x)}, 

2. "vortex-only" configurations, {Z^(x) = sign Tr[U^(x)]}, and 

(R) + 

3. "vortex-removed" configurations, {U^ (x) =Z 1 ll {x)U ll {x)}. 

The procedure of vortex removal, first introduced by de Forcrand and D'Elia [12], is 
known to remove the string tension, eliminate chiral symmetry breaking, and send the 
topological charge of lattice configurations to zero. 

Each of the three ensembles was brought to Coulomb gauge by maximizing, on each 
time- slice, 

^w/(0 = EEiTr[t/,(x,0]. (20) 

x k=\ 




Our results for full configurations at /3 = 2. 1 are exemplified in Figure 1, for a variety 
of lattice volumes. One can observe a sharp "bend" near A — > 0. Both quantities behave 
near like X p , X q , with p, g small numbers. A scaling analysis along the lines of random 
matrix theory (see Appendix B in [7] for details and references) gives the estimates 

A 038 , 



p(A) ~ A 



0.25 



F(X) 



(21) 



with a subjective error of about 20% in the exponents. The Coulomb gauge confinement 
condition, Eq. (16), is thus clearly fulfilled. This result provides a confirmation of the 
mechanism envisaged in the Gribov-horizon scenario. 

To pinpoint the mechanism which might be responsible for the enhancement of 
eigenvalues near the horizon, we considered the Faddeev-Popov observables in the 
vortex-only configurations. Our data are displayed in Figure 2. The enhancement near 
is even more pronounced; it appears that both quantities converge to a non-zero value in 
the infinite volume limit 



p(0)~0.06, F(0)~1.0, 



(22) 



which is confirmed also by a RM scaling analysis. Once again, the confinement criterion 
(16) is obviously satisfied. 



Eigenvalue density, p=2.1 



F(X) s <( Pl l-V 2 | 9x >, p=2.1 
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FIGURE 2. The density p (A) and F(X) for vortex-only configurations. 
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FIGURE 3. The density p (X ) and F (X ) for vortex-removed configurations. 



Finally, results for "vortex-removed" configurations are shown in Figure 3, for the 
largest available, 20 4 lattice volume. The eigenvalue spectrum is strikingly different in 
this case. The eigenvalue density consists of a series of distinct peaks, while values of 
F(X) are organized into separate bands, the lowest few of them clearly separated by 
gaps. This result can easily be understood by inspection of the eigenvalue spectrum of 
the Laplacian operator (equal to M at zero-th order in the coupling). The eigenvalue 
density, at finite volume, is a sum of delta-functions, and each eigenvalue is multiply 
(A^-times) degenerate. The quantity F(Ajt) is equal to Xk- Now if one compares the 
degeneracy of the fc-th eigenvalue with the number of eigenvalues inside the fc-th 
"band" of F{X) (the right panel of Figure 3), one finds a precise match. This can be 
simply interpreted: the vortex-removed configuration is just a small perturbation of the 
zero-field limit 1/^ = 1. This perturbation lifts the degeneracy of degenerate eigenvalues 
and spreads them into bands of finite width. The estimate of the Coulomb self-energy in 
this case indicates that it remains IR finite in the infinite volume limit. 

In pure gauge theory, we used the procedure of de Forcrand and D'Elia [12] to 
study effects of vortex removal on the density of Faddeev-Popov eigenvalues in the 
small X region. However, in some lattice models, confining vortex configurations can be 
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FIGURE 4. The density p (X) for a gauge-Higgs system in the "confined" phase (left) and the "Higgs" 
phase (right). 



suppressed by changing the coupling constants. A prototype example is the gauge field 
coupled to a scalar field of unit modulus in the fundamental representation of the gauge 
group [13]. The action of the model is 

S = S W + |^Tr[0t( x )[/ M ( x )0( x + J u)], (23) 

where Sw is the usual Wilson gauge action and § is an SU(2) group-valued field. In the 
strict sense, this theory is non-confining for all values of couplings /3 , 7, due to the well- 
known Osterwalder-Seiler-Fradkin-Shenker theorem [14, 15]. Even though there is no 
thermodynamic phase transition separating the pseudo-Higgs phase from the pseudo- 
confinement phase of the model, there is a symmetry breaking transition between the 
phases (a "Kertesz" line), and center vortices percolate in the "confinement" phase, and 
cease to percolate in the "Higgs" phase [6, 16, 17, 18]. By varying 7 at fixed gauge 
coupling /3 , one can modify the vortex content and study its effects on Faddeev-Popov 
eigenvalues. 

In Figure 4 I show the eigenvalue density at /3 = 2. 1 for two values of gauge-Higgs 
coupling: 7= 0.6, deep in the "confinement" phase, and 7= 1.2 corresponding to the 
"Higgs" phase. In the former, the densities for full, vortex-only, and vortex-removed 
configurations clearly resemble those in pure gauge theory, while the density in the 
Higgs phase for full configurations looks almost identical to the vortex-removed data 
in the confinement phase. 



COULOMB ENERGY 

We have seen that the Coulomb self-energy of a color non-singlet state is IR divergent, 
due to the enhanced density of Faddeev-Popov eigenvalues near zero. Another question 
is whether the color Coulomb potential of a charge-anticharge pair grows linearly and, if 
so, whether it is also sensitive to the presence/absence of center vortices. This question 
was addressed in Ref. [5], and I will summarize here briefly the main results. 

Let l^qq) = q a (0)q a (R)\ x ¥o) denote a physical heavy static quark-antiquark state in 
Coulomb gauge. Then 

<? qq =(y qq \H\y qq }-(y \H\y }=E se + V coul (R) (24) 

is a sum of self-energy contributions and the 7?-dependent color Coulomb potential. It 
can be computed from the correlator of two timelike Wilson lines in Coulomb gauge: 



G(R,T)=aTr[lS(x,T)L(j,T)]) , ■ L(x,7/)=exp' 



*=l*-y| 



i { dtAo(x,t) 
Jo 



(25) 



It is easy to show that 



£ = E se + V coul (R) = \im V(R,T); S min = E' se + V(R) = lim V(R,T), (26) 
where 

V(R,T) = -£\og[G(R : T)]. (27) 



String tensions extracted from V(R,T) 
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FIGURE 5. cr(r) vs. T . Solid lines indicate the accepted values of the asymptotic string tension at each 
j3 value, with dashed lines indicating the error bars. 



Above, (g'min is the minimal energy of a state containing two static charges, and V(R) 
is the static interquark potential. Since £ > S m i n , it is clear that if V(R) is confining, 
then so is V C0U i(R), and V(R) is bounded from above by V cou i(R), as first proven by 
Zwanziger [19]. 

On a lattice we introduce 



L(x, T) = U (x,a)U (x, 2a)... U (x, T), 



(28) 



and 



V(R,T) = -±log 



G(R,T + a) 



V(R,0) = - 1 -log[G(R,a)}, (29) 



G(R,T) 

from which we obtain an estimate of V cou i(R) (exact in the continuum limit). 

Our results, for SU(2), are shown in Figures 5 and 6. 2 The former figure represents 
a consistency check of our procedure. It verifies that the string tension o(T) extracted 
form V(R, T) approaches the accepted value of the asymptotic string tension at large T. 
In Figure 6 we plot the potential V (R, T) for T = (our estimate of the color Coulomb 
potential) and for T = 4 (approaching the static quark-antiquark potential V (/?)). Upper 
lines of data points in both figures clearly demonstrate that both V C0U i(R) and V(R) 
rise linearly with distance. However, the slope of this linear rise is larger in the color 
Coulomb potential, a cou i 3<7. 3 The fact that the color Coulomb potential overconfines 
does not contradict Zwanziger's bound [19], nor is really surprising. There is no reason 



2 First preliminary results of the determination of the color Coulomb potential in SU(3) lattice gauge 
theory were presented by Nakamura [20] at this conference. 

3 The color Coulomb potential was measured using different methods by Cucchieri and Zwanziger [21], 
and Langfeld and Moyaerts [22], with lower values for the ratio O cou \ja. The origin of this discrepancy 
has not been clarified yet. 




FIGURE 6. V(R,0) and V(R,4) at j8 = 2.5. The former tends to the color Coulomb potential in the 
continuum limit, while the latter approximates the static asymptotic interquark potential. 



to believe that the quark-antiquark state in the Coulomb gauge is the true QCD flux tube 
state, the minimal energy state containing a static quark and antiquark. The Coulombic 
force can be lowered to the true asymptotic one e.g. by constituent gluons present in the 
QCD flux tube, as was suggested in the gluon chain model of Greensite and Thorn [23]. 

The effect of vortex removal is illustrated in lower lines of data points of Figure 6. 
Removing center vortices also removes the confining property of the color Coulomb 
potential. Since the potential, in Coulomb gauge, is sensitive to the interaction kernel K 
in the Hamiltonian (1), and the kernel in turn to the density of near-zero eigenvalues of 
the Faddeev-Popov operator, this result again confirms the observation of the previous 
section: Removal of center vortices alters the density quite drastically. 

We also note that the confining property of the color Coulomb potential is tied to 
the unbroken realization of a remnant global gauge symmetry in Coulomb gauge. This 
connection was studied in detail in Ref . [6] . It was demonstrated there on a few exam- 
ples (deconfined phase in pure gauge theory, pseudo-confinement phase of the gauge- 
fundamental-Higgs theory) that confinement in the color Coulomb potential is not iden- 
tical to confinement in the static interquark potential, and center symmetry breaking, 
spontaneous or explicit, does not necessarily imply remnant symmetry breaking. 



CONCLUSIONS 

Results of our numerical simulations suggest an appealing picture: The low-lying eigen- 
values of the Faddeev-Popov operator in Coulomb gauge tend towards zero as the lattice 
volume increases. The density of the eigenvalues goes as a small power of X, and this, 
together with a similar behavior of the average Laplacian, F(X), assures the infrared 
divergence of the energy of an unscreened color charge. Also, due to the enhancement 
of near-zero modes of the Faddeev-Popov operator, the Coulomb energy of a pair of 
color charges rises linearly with their separation. Both facts support the ideas of the 
Gribov-horizon confinement scenario. 

The constant density of low-lying eigenvalues can be attributed to the vortex compo- 
nent of gauge-field configurations. A thermalized configuration in a pure gauge theory 



factors into a confining piece (the vortex-only part), and a piece which closely resem- 
bles the lattice of a gauge-Higgs theory in the Higgs phase (the vortex-removed con- 
figuration). This establishes a firm connection between the center-vortex picture and the 
Gribov-horizon scenario. This connection is exemplified also by the fact that vortex re- 
moval removes the color Coulomb string tension of the color Coulomb potential. It is 
also consistent with recent investigations of Gattnar et al. [24] in Landau gauge. 

In this talk, I covered results of our numerical investigations. Related analytical de- 
velopments were omitted: (i) thin center vortices lie on the Gribov horizon; (ii) the Gri- 
bov horizon is a convex manifold in lattice configuration space and thin center vortices 
are conical singularities on that manifold; (iii) the Coulomb gauge is an attractive fixed 
point of a more general gauge condition, interpolating between the Coulomb and Landau 
gauges. Interested readers are invited to find them in our recent publications [6, 7]. 
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